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ized deeply inelastic scattering from their Melhn moments to 3-loop order. 



XII Advanced Computing and Analysis Techniques in Physics Research 
November 3-72008 
Erice, Italy 



* Speaker. 

^This project was supported in part by German DFG SFB TR-9, Studienstiftung des Deutschen Volkes, projects 
P19462-N18, P20162-N18, and P20347-N18 of the Austrian FWF, and EC MRTN HEPTOOLS, contract MRTN-CT- 
2006-035505. 



(c) Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



From Moments to Functions 



Johannes Bliimlein 



1. Introduction 

Higher order calculations in Quantum Field Theories easily become tedious due to the larger num- 
ber of terms emerging and the sophisticated form of the contributing Feynman parameter integrals. 
This appUes already to zero scale and single scale quantities. Even more this is the case for prob- 
lems containing more than one scale. While in the latter case the mathematical structure of the 
solution of the Feynman integrals is widely unknown, it is explored to a certain extent for zero- 
and single scale quantities. Zero scale quantities emerge as the expansion coefficients of the run- 
ning couplings and masses, as fixed moments of splitting functions, etc. They can be expressed 
by rational numbers and certain special numbers as multiple zeta-values (MZVs) [1,2] and related 
quantities. 

Single scale quantities depend on a scale z which may be given as a ratio of Lorentz invariants 
s'/s in the respective physical problem. One may perform a Mellin transform over z 



All subsequent calculations are then carried out in Mellin space and one assumes N ^"N, N > 0. 
By this transformation the problem at hand becomes discrete. One may seek a description in terms 
of difference equations. Zero scale problems are obtained from single scale problems treating as 
a fixed integer or considering the Umit N ^°o. 

A main question concerning zero scale quantities is: Do the corresponding Feyimian integrals 
always lead to MZVs? In the lower orders this is the case. However, starting at some order, even 
for single-mass problems, other special numbers will occur [3]. This makes it difficult to use 
methods like PSLQ [4] to determine the analytic structure of the corresponding terms even if one 
may calculate them numerically at high enough precision since one has to known the respective 
basis completely. ^ 

Zero scale problems are much easier to calculate than single scale problems. In some analogy 
to the determination of the analytic structure in zero scale problems through integer relations over 
a known basis (P S LQ) one may think of an automated reconstruction of the all-A^ relation out of a 
finite number of MelUn moments given in analytic form. This is possible for recurrent quantities. 
At least up to 3-loop order, presumably even to higher orders, single scale quantities belong to this 
class. Here we report on a general algorithm for this purpose, which we appUed to the problem 
being currently the most sophisticated one: the determination of the anomalous dimensions and 
Wilson coefficients to 3-loop order for unpolarized deeply-inelastic scattering [8]. Details of our 
calculation are given in Ref. [9]. 

2. Single Scale Feynman Integrals as Recurrent Quantities 

For a large variety of massless problems single scale Feynman integrals can be represented as 

'in a recent analysis [2] the relations between all MZVs in the non-alternating and alternating case were determined 
up to weight w = 12 and up to w = 24 in the non-alternating case using shuffle-, stuffle- and generalized doubling 
relations. These relations lead to basis lengths according to the conjectures by Broadhurst [5] and Broadhurst-Kreimer- 
Zagier [6], which meets the upper bound set by Terasoma, Goncharov and Deligne [7]. In the non-alternating case we 
verified that the latter bound is vaUd at least to vf = 26 and stopped the calculation afterwards due to the large complexity 
involved. All this does not exclude the existence of exotic relations reducing the basis further. 
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polynomials in the ring formed of the nested harmonic sums Sai,...,a,,{N), [10, 11] and the MZVs 
over the rational function field Q{N). Here, 



^W(A^) = I ^^^^P^S.ik), ai,b G Z\{0}, (2.1) 
Cm = I ^^^^S,ik), aub G Z\{0} . 

If ft = 1, the meaning of i^^^ is symbolic, since it diverges. The degree of divergence is a positive 
power of S\ (oo). Rational functions in A'^ and harmonic sums obey recurrence relations. Thus, due 
to closure properties [14] also any polynomial expression in such terms is a solution of a recurrence. 
Consider as an example the recursion 



(A^+l)l"l 

It is solved by Sa{N). Corresponding difference equations hold for harmonic sums of deeper nest- 
edness. Feynman integrals can often be decomposed into a combination containing terms of the 
form 

dz-^ H^aiz), / ^Z^4t ^a{z) , 

1— Z JO 1+z 

with H^{z) being a harmonic polylogarithm, [12]. This structure also leads to recurrences, cf. [13]. 
It is very hkely that single scale Feynman diagrams do always obey difference equations. 

3. Establishing and Solving Recurrences 

We assume that a sufficiently large set of moments at integer values A'^,- is given for a physical 
quantity, which obeys a recurrence relation. One seeks 



k=Q 



d 

[ 

1=0 



F{N + k)=0. (3.1) 



The method for determining potential recurrences is available in standard packages [14]. The cor- 
responding linear system is dense. Rational number arithmetics is not feasible for the large systems 
to be solved. Let us consider as an example the difference equation being associated to the contri- 
bution of the color factor Cp for the 3 -loop Wilson coefficient ^ in unpolarized deeply inelastic 



scattering. 1 1 Tb of memory would be required to establish ( p.l\ ) in a naive way. Therefore refined 
methods have to be applied. We use arithmetic in finite fields together with Chinese remainder- 
ing [15], which reduces the storage requirements to a few Gb of memory. The linear system 
approximately minimizes for I ^ d.lf one finds more than one recurrence the different recurrences 
are joined to reduce / to a minimal value. It seems to be a general phenomenon that the recurrence 
of minimal order is this with the smallest integer coefficients, cf. also [16]. For even larger prob- 
lems than those dealt with in the present paper, a series of further technical improvements may be 
carried out, [17]. 
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For the solution of the recurrence low orders are clearly preferred. It is solved in depth-optimal 
nZ fields [18]; here we apply advanced symbolic summation methods as: efficient recurrence 
solvers and refined telescoping algorithms. They are available in the summation package Sigma 
[19] implemented in the computer algebra system Mathematica. 

The solutions are found as linear combinations of rational terms in N combined with functions, 
which cannot be further reduced in the nZ fields. In the present application they turn out all to be 
nested harmonic sums S-^{N), (2J_). Other or higher order applications may lead to sums of different 
type as well, which are uniquely found by the present algorithm. 



4. Determination of the 3-Loop Anomalous Dimensions and Wilson Coefficients 

We apply the method to determine the unpolarized anomalous dimensions and Wilson coefficients 
to 3-loop order. Here we use the above method to all the contributions to a single color/i^, -factor. 
These are 186 terms. As input to the calculation we use the respective Mellin moments, which 
were calculated by a MAPLE-code based on the harmonic sum representation [8]. We need very 
high moments and calculate the input recursively. As an example, let us illustrate the size of the 
moments for the C| -contribution to the Wilson coefficient C2.q{N). The highest moment required is 
N = 5114. It cannot be calculated simply with summer [1 1] because of time and storage reasons. 
The highest moment is a rational number with a numerator of 13388 and a denominator 13381 
digits. Below we give the moments for A/^ = 3 and N = 500. 

N-3 : 

#11 digits / #10 digits 
-98268084191 / 1166400000 
N-500 : 

#1262 digits / #1256 digits 

16418407704241967 80 953020 61917 637 650 62843035444812 62083057197 6007465070084 937 93 994 
42241103234415 91630311482222058287 688 94220 957085 9151121677307585313995100 9783 63179 
2518952817622034037186132846974627021672678012913675099511203807811938593043910803 
5044345 92021869605258833203 63553250899983 6135422 68823 6732214 90376310537617 64348772 
54038108742 64 96872 9520075 6192272854718024194037272078224737 65 9999002 36383740315299 
2050533 60163348434824 94547575553446642108141111400 6547539113 67 98 6891674100 6507 6749 
357870 9478 6833 90573 9774100135208 944 944 63 90 92 913274258157665 66386397276158317387748 
5 9454713 92 64 608 97008751574450750731923285428 90 9654 62004805711998 74814441437 938 60 93 
99373 617 9802 904442578995372 6133 6751997 905237704272 985005100 634 64 0619858400662 96071 
3372543015 64 8 91915596406960 69945 973638863011850 678272 919370 653007547869470 63 672848 
938208192 6871078 600328 62813193 67 66057475 970450896556667 6221633 65895808773428119721 
53527 921310890 635770450 696939622130 611988 94057033 60 60 68 695 6071232719697269810 60056 
011584 60943 6023 998 62339178722 60722277322 69045013237 683 625354 9152130116645 6705 65045 
9666945 92016458 6023 9580 6027174660 67 988 988 613 607723330880307417756055465187887 93327 
22 643 682 97071217405 65447437584423825088 9238538 9745484212 9817042590952174255 94 94728 
72017877003 9473965 622616598603668391544078534 62338171648227013134266795320251847 
/ 

3057444 6142472253728 825705143 6735 8 6 9727 813074134 82 8212220 64 92 932 82 0352 4 40 850471902 
74 9104 696210533 66455 63 654873 675 6907 96713 90 65 65 6888203 65 6019072 6371086395482 638 6081 
322758003787 93 618 699410038028075908 603588 941428 9104 677 6447162895908787986423254 678 
577 6778283337231702130 61249942 981955 97 98501074020 67 628276928 91029556794218857958 67 
1982932 998 601320344 971927374 90588 993405 998727193 97 6021283 63 68 61950118 92382154423 66 
380577370192 950 92 6815774799285 93848374037511830194236928 68569168206789710047557452 
51312173822720 60267681480496298 9755224 67 61470784863 9773185 90 9858278799786 637303834 
101716667 627 68475257047554 931662 632 9707 9720470719813 623 90154581195385398 645 6533543 
99941820505518279599887 6012116849074547 696925 94 68454 61343162417 91988607 5151307 6481 
030473420592 6703138519418575731315944374873897873 64 670 6993 620825697218523316375559 
40 682220047 65 962715 92420852 610 600810974040238012 62 60947524 640509361283802755722132 
4856690051525724 685 9197 92 64150 60823075 67 95 69623285 6007347108 67 992 871312 875 64 668441 
62566980835042338 9743 64847020024713143308034214 67773 925541151273 92 498594 6178771189 
23124371622134381377038960 6473498715702080141315355543531132 6719739117599044341913 
5 922693587373856 60 95 94245 94 823746 92 931487025167140382 97077 63 93 82 3322512 55360181047 
4965862324750 911265 97 62997 67 973752788271111167745 9300352000000000000000000 

N-5114 : 

#13388 digits / #13381 digits 
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The corresponding difference equations ( ^ ) are determined by a recurrence finder. Furthermore, 
the order of the difference equation is reduced to the smallest value possible. The difference equa- 
tions are then solved order by order using the summation package Sigma. 

For the C|--term in C2,q{N) the recurrence was established after 20.7 days of CPU time. Here 
4h were required for the modular prediction of the dimension of the system, 5.8 days were spent 
on solving modular linear systems, and 11 days for the modular operator CCDs. The Chinese 
remainder method and rational reconstruction took 3.8 days. 140 word size primes were needed. As 
output one obtains a recurrence of 3 1 Mb, which is of order 35 and degree 938, with a largest integer 
of 1227 digits. The recurrence was solved by S igma after 5.9 days. We reached a compactification 
from 289 harmonic sums needed in [8] to 58 harmonic sums, where the representation in [8] (see 
a corresponding attachment) has already been compactified following an idea of one of the present 
authors. The determination of the 3-loop anomalous dimensions is a much smaller problem. Here 
the computation takes about 1 8 h only for the complete result. 

For the three most complicated cases, establishing and solving of the difference equations took 
3 + 1 weeks each, requiring < 10Gb on a 2 GHz processor. This led to an overall computation time 
of about sixteen weeks, with the possibility to parallelize four times. Here we did not yet consider 
parallelization w.r.t. the 140 primes chosen, which would significantly reduce the computational 
time of the C| term discussed above and for other comparably large contributions. 

In the final representation, we account for algebraic reduction [20]; for this task we used the 
package HarmonicSums [21] which complements the functionalities of Sigma. One observes 
that different color factor contributions lead to the same, or nearly the same, amount of sums at 
a given quantity. This points to the fact that the amount of sums contributing, after the algebraic 
reduction has been carried out, is governed by topology rather than the field- and color structures 
being involved. The linear harmonic sum representations used in [8] require many more sums than 
in the representation reached by the present analysis. A further reduction can be obtained using the 
structural relations, which leads to maximally 35 different sums up to the level the 3-loop Wilson 
coefficients [13]. It is not unlikely that the present method can be applied to single scale problems 
in even higher order As has been found before in [13,22-24] representing a large number of 2- and 
3-loop processes in terms of harmonic sums, the basis elements emerging are always the same. This 
applies to the anomalous dimensions and Wilson coefficients of the space- and time-like polarized 
and unpolarized case, the polarized and unpolarized Drell-Yan process and hadronic Higgs-Boson 
production in the heavy mass limit, deep-inelastic heavy flavor production in the limit ^ nig, 
higher order QED corrections in e^e^ annihilation, as well as to soft and virtual corrections to 
Bhabha scattering. 

In practice no method does yet exist to calculate such a high number of moments ab initio as 
required for the determination of the all N formulae in the 3-loop case. On the other hand, a proof of 
existence has been delivered of a quite general and powerful automatic difference-equation solver, 
standing rather demanding tests, which is ready to help in the solution of present day problems in 
higher order Quantum Field Theory. It opens up good prospects for the development of even more 
powerful methods. 
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5. Conclusions 

We established a general algorithm to calculate the exact expression for single scale quantities from 
a finite (suitably large) number of moments, which are zero scale quantities. The latter ones are 
much more easily calculable than single scale quantities. We applied the method to the anomalous 
dimensions and Wilson coefficients up to 3-loop order. To solve 3-loop problems this way is not 
possible at present, since the number of required moments is too large for the methods available. 
We established and solved the recurrences for all color resp. (^-projections at once, which forms 
a rather voluminous problem. Yet we showed that giant difference equations [order 35; degree 
~ 1000] can be reliably and fast established and solved unconditionally for the most advanced 
problems in Quantum Field Theory. 

Acknowledgment. We would like to thank J. Vermaseren for discussions. 
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